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Abstract. We describe the results of a search for transit-like events caused by giant planets occulting stars in the 
old, metal-rich open cluster NGC 6791 based on BV^J-colour photometry from eight nights of observations with 
the 2.54 m Nordic Optical Telescope. To extract the light curves we have used both PSF photometry (daophot) 
and the difference imaging technique (isis). We have re-analyzed observations from earlier campaigns to search 
for multiple transits, determine periods of long-period variables, and detect eclipsing binaries. We confirm 20 
known variables and have discovered 22 new low-amplitude variables with amplitudes in the range 7-40 mmag 
and periods 1-16 days. We have found the primary and secondary eclipses of two eclipsing binaries based on 
the new and older photometric campaigns. The search for transits-like events has turned up a few single-transit 
candidates. The transit depths are 10 mmag in both the V and / filters over periods of 1.0-2.5 hours, but future 
observations are required to see whether identical transit-like events recur in these same stars. 
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1. Introduction 

Since the discovery of the first extra-Solar giant planet 
around 51 Peg (Mayor & Queloz [T995|l . around 100 extra- 
Solar planets have been discovered. Nearly all of these 
have been found by measuring the variation in radial ve- 
locity caused by the orbit of the planet around the parent 
star. In the case of HD 209458, the planet is known to tran- 
sit the star causing a reduction of 1.7% in the light from 
the star. From the photometric and radial velocity curves 
Charbonneau et al. (2000) could determine the mass and 
radius of the planet to within a few percent. Interestingly, 
the resulting mean density of the planet is quite low due 
the close vicinity of the parent star (Burrows et al. ..2000J . 

Several attempts have been made to detect extra-Solar 
planets by monitoring many thousand stars simultane- 
ously with an imaging CCD camera in order to detect the 
reduction in light due to the transit of a planetary com- 
panion. Most notable is the result for 47 Tuc (Gilliland 
et al. 2000), which failed to detect planets during an 
8.3 day search with the Hubble Space Telescope. 



Send offprint requests to: H. Bruntt, bruntt@phys.au.dk 

* Based on observations made with the Nordic Optical 
Telescope, operated on the island of La Palma jointly by 
Denmark, Finland, Iceland, Norway, and Sweden, in the 
Spanish Observatorio del Roque de los Muchachos of the 
Institute de Astrofisica de Canarias. 



In addition, Udalski et al. (2002a) selected a sample of 
52 000 stars from the OGLE-III (Wozniak et al. 2002) cam- 
paign. These stars have photometry better than 1.5% and 
were observed over 32 nights during a period of 45 days. 
Udalski et al. reported 47 transiting systems of which 
42 exhibited multiple transits. An additional 13 objects 
were later reported by Udalski et al. (2002b), who used 
a refined search method. In a subsequent paper Dreizler 
et al. (2002) presented spectra of 16 of the most promising 
planet transit candidates found by Udalski et al. (2002a). 
From these spectra Dreizler et al. (2002) could determine 
the spectral classes of the stars and in turn constrain the 
mass and radius ratios of the orbiting object and the par- 
ent star. Dreizler et al. found that only two of the 16 tran- 
siting objects have radii and masses consistent with being 
giant planets. The other transiting objects are probably 
brown dwarfs or M-type stars. 

Gonzalez et al. (2000) have made spectroscopic stud- 
ies to determine the metallicities of stars with planets and 
from a sample of 39 stars they find <[Fe/H]>= 0.18±0.19. 
In comparison, the local average metallicity of G-type 
dwarfs is [Fe/H] = -0.2 dex (Rocha-Pinto & Maciel 1996). 
Laughlin (2002) states that the probability of a star host- 
ing a planet is 10 times higher for stars with high metal- 
licity ([Fc/H] > 0.2) than for less metal-rich objects. One 
possible reason for the higher metallicity of stars with 
planets could be that planets are more likely to form in 
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a metal-rich environment (see Laughlin 2002 for a discus- 
sion). Another possibility (Laughlin 2002) could be the 
addition of metal-rich material onto the star, which is 
mixed up in the convective outer layers. The source could 
be migrating proto-plancts and/or planetesimals from the 
proto-planetary disc. 

Laughlin (2000) also examined the correlation between 
metallicity and mass of stars with giant planets. Among 
stars with planets there is a tendency for the stars with 
higher metallicity to also have larger mass. The proposed 
explanation is the fact that the higher the mass the shal- 
lower the outer convection zone will be, hence for the more 
massive stars the metal-rich material which is added is 
mixed through a smaller fraction of the star. 

Inspired by these interesting results we selected the 
old metal-rich open cluster NGC 6791 as our target for 
a search for planet transits. The metallicity of the clus- 
ter is a subject of debate. Taylor (2001) states that "the 
best that can be said at present is that the cluster metal- 
licity lies in the range from -1-0.16 to -fO.44 dex". Two 
studies published since the Taylor article come to no bet- 
ter agreement: Friel et al. (2002) have derived [Fe/H] = 
-1-0.11 ± 0.010 from spectroscopy of cluster giants, while 
Stetson et al. H2003|l infer an abundance ^ +0.3 by match- 
ing the morphology of the cluster color-magnitude dia- 
gram to that of recent theoretical isochrones. 

In comparison the globular cluster 47 Tuc studied by 
Gilliland et al. (1999) has [Fe/H] = -0.7 (Salaris & Weiss 
1998). In this study around 30 transits of short-period gi- 
ants were expected from the 34000 light curves, i.e. one 
out of every 1 000 stars. The fact that no transits were 
found is perhaps due to the lower probability of finding 
planets in a metal-poor system. Another possible expla- 
nation is that the relatively high stellar density in a glob- 
ular cluster may have influenced the evolution of proto- 
planetary discs or the long-term stability of any planetary 
systems that may have formed. The dynamical properties 
of open clusters are not as extreme as for globular clusters 
and the formation of planets may take place here. 

We observed NGC 6791 for ~ 8 hours on each of seven 
nights and we estimate the probability of detecting a sin- 
gle transit of a close-in planet to be Psingic — 50% and 
Ptwice — 7% for detecting a transit twice for the same star. 
The probability that a local G or F-type field star har- 
bours a close-in giant planet is around 0.7 % (Butler et al. 
I2002[l while the geometric probability that the planet will 
cause a transit is 10%. We have -/Vct<2% = 2 500 stars in 
our sample with light curves with a precision better than 
2% per data point. Assuming that the binary fraction Pbin 
of the cluster is 15% the expected number of planets that 
will be detected is 

N = iV^<i% * PgF * Pgoom * (1 - Pbin) * Psinglo = 0.8 (1) 

The metallicity of NGC 6791 is significantly higher 
than the field stars and we may expect to find at least a 
few candidates with single transits. According to Laughlin 
H2000|l the fraction of stars harbouring planets is ~ 10 
times higher for metal rich stars. 



If we detect transits in NGC 6791 and combine this 
with radial velocity measurements we can put constraints 
on the mass and radius of the parent star and in turn 
also constrain the radius of the transiting object from the 
depth of the transit. 

In Sect. 2 we describe the observations and in Sect. 3 
we discuss the data reduction. In Sect. 4 we present the 
calibrated colour-magnitude diagram. In Sect. 5 we discuss 
the detected transit candidates. In Sect. 6 and 7 we present 
the new variable stars and the detached eclipsing binaries. 
Finally, in Sect. 8 we discuss our results and suggest future 
prospects for planet searches though photometry. 

2. Observations 

The observations for this project were acquired between 
July 9 and 17, 2001 at the Nordic Optical Telescope 
(NOT) by one of us (Bruntt). We used the ALFDSC instru- 
ment, which gives a field of view of 6' 5 x 6' 5. A 2k x 2k 
thinned Loral CCD with a pixel size of 0"188, read-out 
noise 6.5 e~ (rms), and a gain of 2.5 e~/ADU was used 
as detector. The time-series observations were made us- 
ing V and / filters to be able to distinguish between grey 
and coloured transits. To first order, planets will cause the 
same transit depth regardless of colour, while a transit due 
to a stellar companion or variations due to the star itself 
(spots) will most likely have different amplitudes in dif- 
ferent filters. The typical exposure times used were 180 s 
and 90s for V and / respectively, chosen such that stars 
above the turnoff of the cluster at around V = 16 would 
be at the saturation limit. In order to reduce the overhead 
in connection with changing filters we grouped the obser- 
vations as five exposures in / followed by three in V. The 
read-out time of the detector was 70 s. 

A few observations were collected in the B filter as well, 
to be used for the construction of the colour-magnitude 
diagram of the cluster. 

In order to achieve the highest possible precision in the 
differential photometry we made small corrective offsets 
of the telescope during the night to keep the stars in the 
same position on the detector. To be specific, the l-cr RMS 
scatter of the x and y positions was 2.8 pixels and the 
median seeing was 5.3 pixels. 

Flat fields were obtained with different rotations of the 
camera (to minimize scattered light effects; see Grundahl 
& S0rensen I1996|l during evening and morning twilight 
on each night. A master fiat field was constructed by us- 
ing all flat field images, since no significant night-to-night 
variations were seen. 

On each night we observed several standard fields from 
Stetson l|2000|l in order to provide a good transformation 
to the standard BVI system^ (Johnson BV and Cousins 
/). The details of this can be found in Stetson et al. H2003fl . 
For the remainder of this paper we will refer to the data 
set from this run as NOTOl. 

^ Available from |http : / /cadcwww . hia . nrc . ca/ cadcbin/ 1 
wdb/ astrocat/ stetson/ query/ . 



H. Bruntt et al.: A Search for Planets in the Old Open Cluster NGC 6791 



3 



Table 1. Fundamental properties of NGC 6791: right ascension, declination, age, distance modulus, metallicity, and 
reddening. The values are taken from Chaboyer et al. (,1999J except the metallicity which is from Taylor ( 20O1J • Note 
that the quoted errors are probably too optimistic and do not include systematic errors (except the error on the 
metallicity) . 

a2ooo.o <52ooo.o Age [Gyr] {m ~ M)v [Fe/H] E{B-V) 

19''20™53" +37° 46 30 8.0(5) 13.37(10) +0.30(15) 0.10(2) 



The fundamental parameters of NGC 6791 are given 
in Tabled 

2.1. Time-Series Studies of NGC 6791 

Rucinski et al. H1996|l carried out a campaign on 
NGC 6791 for 14 nights in July 1995 to search for de- 
tached eclipsing binaries. Despite good weather conditions 
the noise level in this data set (labelled SMR95 from this 
point on) is higher compared to the NOTOl data due to 
the smaller telescope: the apertures are 1 m for SMR95 
and 2.5m for NOTOl. The SMR95 data is still useful 
for constraining the periods of long-period variables and 
will be used to search for detached eclipsing binary stars. 
Rucinski et al. H1996|l discovered 11 variable stars (four 
new ones) and among them were three detached binaries. 
We refer to Rucinski et al. H1996fl for a list of earlier stud- 
ies on NGC 6791. 

A more recent study was made by Mochejska et al. 
H2UU2|I who have selected NGC 6791 as their first target 
to search for planetary transits. They have used the FLWD 
1.2 m telescope at Mt. Hopkins, Arizona, from July 6 to 
August 1, 2001 (we label this data set MOCHOl), i.e., 
partly overlapping with the NOTOl observations from La 
Palma. Unfortunately, the weather conditions were quite 
poor with a median seeing in V of 2''1 (in contrast, for 
NOTOl the median seeing was 1''0). Mochejska et al. 
H20()2|l state that the search for transit-like events still re- 
mains to be done, and that a new campaign on the cluster 
is being planned. Mochejska et al. I)2()[)8|l have analysed 
time-series data spanning six years and report the dis- 
covery of seven new long-period and irregular variables. 
Kaluzny |2003l reported the discovery of four new vari- 
ables in the field. 

In Table |21 we give a summary of the observations of 
NGC 6791 that we have used. For the MOCHOl data 
we have only used the light curves of the variable stars 
from Mochejska et al. H2()()2|l : in other words, we have 
not reduced these data ourselves as we have for NOTOl 
and SMR95. Note that the photometric precision of the 
NOTOl campaign is superior to the other telescopes, and 
only these data (seven nights) were used to search for 
transit-like events. 

In addition to the data listed in Table|21several images 
taken outside the center of NGC 6791 from NOTOl and 
several images from other telescopes were reduced with 
ALLFRAME as described by Stetson et al. H2008|l . 



3. Photometry of NGC 6791 

We have used two different methods for obtaining the 
differential time series photometry: (1) the profile fit- 
ting (PSF) approach viz. daophot/ allstar/ allframe 
(Stetson 119871 11994(1 and (2) the image-subtraction 
method using the ISIS software (Alard & Lupton 11998(1 . 
We shall describe our approach for these two methods be- 
low. 

3.1. ALLFRAME Reductions 

The time-series photometry from SMR95, the new NOTOl 
data set, as well as a few frames from other campaigns 
were reduced using the DAOPHOT / allframe photomet- 
ric software fStetson ll987[ll9i?l|l . This is described in our 
companion paper (Stetson et al. 12003(1 . The individual 
light curves of the variable stars (and the transit-like can- 
didates found from the ISIS light curves) were derived from 
these data as follows. 

Transformation equations relating the instrumental 
CCD magnitudes to standard-system magnitudes and 
colours, with corrections for atmospheric extinction, were 
derived from standard fields observed the same nights as 
the NGC 6791 observations were obtained. These trans- 
formations must then be used again in the reverse sense 
to convert the observed instrumental magnitudes for the 
cluster stars to magnitudes on the standard system. Since 
these transformations are expressed in terms of standard, 
not instrumental colours, initially each given star was as- 
sumed to have a "typical" colour. Then, robust least- 
squares methodology was used to find the standard-system 
magnitude in each of the photometric band-passes which 
was most consistent with all the observed instrumen- 
tal magnitudes for that star. These new standard-system 
magnitudes then implied a new, more realistic colour for 
the star. These were used for a new least-squares solution 
for the best standard-system magnitudes, and the pro- 
cess was iterated until the derived magnitudes stopped 
changing. Because the colour terms in the transforma- 
tions are quite small, with absolute values of well under 
0.1 mag/mag, this process converges very quickly. This ap- 
proach does have the drawback that at each epoch the 
instrumental magnitude is transformed to a standard- 
system magnitude based upon the mean colour of each 
star, averaged over all epochs. In the case of a constant 
star, this is an optimal approach, but in the case of a vari- 
able star it could introduce small systematic errors in the 
magnitudes observed at the extremes of the star's colour 
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Table 2. Log of observations of NGC 6791 from three campaigns. The columns contain the adopted abbreviation of the 
campaign, the name of the telescope, the duration of the campaign, the number of useful nights (n), the median seeing 
in the V filter, and the number of images obtained in each of the BVRI filters. The last column gives information on 
which reduction technique we have used. Note that we have not reduced the MOCHOl data set, but we have obtained 
the light curves of the variable stars published by Mochejska et al. I|2()()2|l . 



Campaign 


Telescope 


Dates UT 


n 


["] 


Number of imaj 
Nb Nv Nr 


^es 
Ni 


Reduction 
Method 


NOT 01 


2.5 m NOT 


9/7-17/7 01 


7 


l'!0 


58 264 


AAA 


allframe/isis 


MOCHOl 


1.2m FLWO 


6/7-1/8 01 


18 


2'^1 


36 204 




Not Reduced 


SMR95 


1.0m JKT 


11/7-24/7 95 


14 


l'^6 


601 


153 


ALLFRAME 



cycle. For a variable star, it would be better to correct each 
instrumental magnitude according to the star's instanta- 
neous colour, which is possible when light curves are be- 
ing fitted at the same time as the photometric observations 
are being calibrated. This was done, for instance, when we 
used the same software for the Hubble Space Telescope Key 
Project on the Extra-galactic Distance Scale, where tem- 
plate Cepheid light curves were being fitted to variable- 
star candidates (see Stetson il9961 [T^M|) . In the present 
case, we had no preconceived notions as to the forms of 
individual variable stars' light curves, so we did not try to 
fit light curves during the photometric reductions in order 
to be able to predict instantaneous colours; the extraction 
of light curves was done ex post facto after the calibration 
of the photometry. Note that this procedure may bias in- 
stantaneous standard-system magnitudes by amounts not 
larger than the colour-transformation terms times half the 
colour amplitude of a variable star, which will generally 
lead to a maximum systematic error less than 0.01 mag in 
the instantaneous magnitudes at the epochs of extreme 
colour. Systematic errors in the periods, mean magni- 
tudes, and mean colours of the variable stars should be 
negligible. 

3.2. Difference Imaging witli isis 

We have used the difference-imaging software ISIS version 
2.1 developed by Alard & Lupton (1998). For each of the 
filters V and /, we selected the images with the best see- 
ing, which were then averaged to generate reference im- 
ages. For each original image a kernel is computed which 
describes the variations of the PSF across that image rel- 
ative to the reference image. The final step is to convolve 
the reference image by this kernel and subsequently sub- 
tract the convolved reference image from each individual 
image. In principle what will be left in the subtracted im- 
ages will be the signal that is intrinsically different from 
the reference image, e.g., variable stars. 

One of the advantages of using the difference-image 
technique is that the signal from variable stars in crowded 
regions will be less affected by the neighbouring stars, e.g., 
when compared to profile fitting photometry. Another im- 
portant point is that variations due to airmass and trans- 
parency variations are removed to first order as a part of 
the image subtraction. 



We tried to use the photometry package of ISIS on 
the subtracted images. Unfortunately, the program pro- 
duced only meaningless output (all stars were constant). 
Instead we used the daophot package to extract aper- 
ture photometry in the subtracted image around the 
position of all stars in the field (found by using ALL- 
STAr). Normally, daophot calculates the magnitude of 
a star as m = 25 — 2.51og(A^), where N is the num- 
ber of counts in the stellar profile. But since the num- 
ber of counts in an aperture can be negative for the dif- 
ference images, we had to add a short piece of code in 
the original daophot program to account for this. In 
particular, we defined a new "difference magnitude" as 
TOisis = 25 - 2.5 log \{N -\- ^N'^ + N^)/2 \ . The zero- 
point of this magnitude scale is mo = 25 — 2.51og(iVo) 
and we used A'o = 1, i.e., m = 25.0 when N = 0. 
One can show that the number of counts will be iV = 
iVosinh[0.41n(10)(mo - m + 2.51og(2))]. 

We used a range of increasing aperture sizes which 
were scaled with the seeing (i.e., the measured fwhm) of 
each image. We then defined magnitude intervals with a 
width of 1 magnitude and selected the aperture size to 
be used within each interval which gave the lowest noise 
level. Large apertures were used for the bright stars and 
smaller apertures for fainter stars (the latter will be more 
affected by sky background noise). Based on growth curves 
constructed for the direct images f Stetson I1990|l we then 
scaled the counts in each aperture to include all the light. 

The relative magnitudes were then computed as Am = 
—2.5 log[(A'rcf — Ni)/N,:ci], where Nrd is the counts in the 
reference image and Ni is the number of counts in the cho- 
sen aperture after applying the correction from the growth 
curve. 



3.3. Repairing the Images 

Any major defects in the images will cause residuals in 
the subtracted images. There are three major contribu- 
tors to the defects: cosmic rays, bad columns, and satu- 
rated stars, including charge leakage extending down the 
image columns from the brightest saturated stars. While 
cosmic rays are normally rejected as part of the ISIS re- 
duction the latter two effects may infiuence the quality of 
the image subtraction. In particular the image subtraction 
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gives large residuals around the most saturated stars. We 
thus decided to repair the saturated pixels in the cores 
of bright stars and the flux from stars which fall on bad 
columns. As a result of this procedure, which we will now 
describe, the residuals in the subtracted images were still 
present, but only very close to the saturated star. 

A short code to replace pixels on the bad columns and 
the saturated pixels was written in idl. For each image 
the code examines the pixels around the brightest stars to 
look for pixels above the saturation limit (60 000 ADU). 
We also identified stars for which the centroid is within 
20 pixels of a bad column. We then computed an "arti- 
ficial" image where profiles for these stars were inserted 
using the PSF model found with daophot for each image, 
employing for this purpose the positions and instrumental 
magnitudes found for these stars by ALLSTAR. We then 
replaced the saturated pixels and the signal on the bad 
columns in the observed image with the data from the 
artificial image (including the estimated background level 
and photometric noise). 

We stress that when searching the light curves for 
transit-like events and variable stars, we clearly mark the 
points in the light curves where the star is close to a bad 
column or a neighbouring saturated star, since the validity 
of the measured signal in these cases is indeed question- 
able. 

3.4. Photometric Noise in tlie isis Light Curves 

In Fig. n] we plot the image-to-image noise level of the 
light curves from the NOTOl data set reduced with ISIS. 
The dashed lines mark the contribution to the noise from 
photon statistics and aperture noise (flat field error, sky 
background determination, and read-out noise) calculated 
using Eq. (31) in Kjeldsen & Frandsen (|1992j): 



2 In 2 
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-I- nr. 



6sky 
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where cr^p the theoretical total noise level, W is the 
FWHM of the stellar profile (in pixels) , eg is the number 
of electrons in the flat field in one pixel, Cstar is the total 
number of electrons from the star, rap is the radius of the 
aperture (in pixels), Csky is the number of electrons per 
pixel in the sky background, and (Tccd is the readout noise 
in electrons. 

The solid line gives the combined theoretical noise 
level. The brightest stars in the cluster reach a noise level 
of Immag per point while the turnoff stars (Jinst = 12) 
typically have 3mmag per point. Of the 8455 stars in our 
sample, 1 700 stars have noise < 1 % and for 2 500 stars 
the noise is < 2% per data point. 

The fact that the measured noise level is much higher 
than the photon noise level for the faint stars shows that 
the combination of image-subtraction and aperture pho- 
tometry on the subtracted images is indeed not yet opti- 
mal. 

The relative transit depth caused by a close-in giant 
planet is equal to the square of the radius ratio of the 
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Fig. 1. Measured noise level in V filter (instrumental mag- 
nitudes) measured over seven nights from the NOTOl data 
set. The reductions were made with ISIS (Alard & Lupton 
.1998 ). The dashed lines mark the photon noise and the 
aperture noise while the solid line is the quadratic combi- 
nation of these two. 

planet and the parent star. For stars below the turnoff 
region (/instr — 12), the transit depth of a Jupiter-sized 
planet will be of the order 1-3% (cf. Fig. 2 in Mochejska 
et al. 2002). Thus the photometric noise limit is well below 
the limit for detection of transits since we will have several 
measurements during a transit which is expected to be of 
the order 1-3 hours. 

4. The Calibrated Colour-Magnitude Diagram 

In Fig. 121 we present the colour-magnitude diagram of 
NGC 6791. It is the result of the ALLFRAME reduction of 
958 V and 737 / frames which is described in our compan- 
ion paper (Stetson et al. I2003|l . The magnitudes are cali- 
brated to the standard V and / system. The formal error 
on the turnoff stars is around 0.5mmag (standard error of 
the mean) in both V and / for stars around the turnoff 
while the calibration errors are probably not much greater 
than 2 mmag. It should be noted, however, that these are 
repeatability estimates for the magnitudes and colours of 
individual stars contained within this data set. Stetson 
et al. present evidence that two stars that appear to have 
identical photometric properties with one filter/detector 
combination can differ by up to 20 mmag (standard devi- 
ation) when measured with different equipment, in these 
broad bandpasses. Nevertheless, the calibrated photome- 
try could be very useful as a standard field and is available 
from the Canadian Astronomy Data Centre (see footnote 
on PageEJ. 

5. The Search for Transit-like Events 

The first seven nights of the NOTOl campaign were used 
in the search for transit-like events. The quality of the 
eighth night was poor by comparison and neglected for 
that reason. Moreover, the data from the SMR95 cam- 
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1.0 1.5 2.0 2.5 

V ~ Ic 

Fig. 2. The calibrated colour- magnitude diagram of 
NGC 6791. The inserted plot shows the turnofF region: 
notice the narrow main sequence and the binary sequence. 



paign was also considered to be too poor for the detection 
of low-amplitude transits. The NOTOl light curves were 
searched for low-amplitude transits using the matched fil- 
ter technique, first suggested for this task by Jenkins et al. 
H1996|l . The decision to use this method to identify transit- 
like features is based on the work of Tingley ;'2()().'i}. where 
it is shown that the matched filter performs best. The data 
from NOTOl was of very good quality and very consistent, 
with no variation of the signal-to-noise ratio from night to 
night. 

Our first implementation of this transit-detection al- 
gorithm searched for single transits without using point- 
to-point noise, which is reasonable considering the consis- 
tency of the data. It assumed a square-well transit shape 
and looked for transits with durations of 2, 3, and 4 hours 
and depths of 1, 2, and 3%. A histogram was made of 
the times of the derived events, to determine if any par- 
ticular epochs were over-producing anomalies, which in 
theory should be distributed evenly in the case of uncor- 
related noise. The epochs that over-produced events were 
then removed and the analysis repeated. This process was 
performed several times, which significantly reduced the 
number of false alarms. Many of the events that remained 
were also easily identifiable as false alarms, as they could 
be explained by the contamination of a star by a nearby 
known variable star, a saturated neighboring star, or by 
the star being close to the edge of the CCD or a bad col- 
umn. However, many interesting low-amplitude variable 
stars were identified, including ones which were not clearly 
identified by the Stetson J-index (cf. Sect. EI). 

The second implementation of the transit-detection al- 
gorithm searched for multiple transits using point-to-point 
noise. The point-to-point noise was calculated by choosing 
a sample of some tens of stars in a relatively narrow magni- 




0.8 0.9 1.0 1.1 1.2 

Y - Ic 

Fig. 3. Colour-magnitude diagram of NGC 6791 with 
standard V and V — Ic magnitudes from Stetson et al. 
(|2003). The transit candidates (triangles) and the eclips- 
ing binaries (circles) are marked. The stars with the best 
photometry are also plotted to emphasize the location of 
the main sequence, the turnoff, and the red giant branch. 



tude bin with similar, low signal-to-noise ratios. These val- 
ues were then weighted by the the overall signal-to-noise 
of the individual light curves. Again using the matched fil- 
ter with a square-well transit shape, the light curves were 
searched for repeating transits with periods ranging from 
0.8 to 6.3 days, durations of 2, 3, 4, and 5 hours, and 
depths of 0.02, 0.01, and 0.005 magnitudes. This imple- 
mentation identified fewer strong events, but fewer false 
alarms and variable stars as well, the latter due to the fact 
that most of them were excluded from the experiment due 
to their relatively high signal-to-noise ratios. 

5.1. Transit-like Light Curves 

After rejection of obvious false alarms, we have found ten 
stars with transit-like events in their light curves. The 
positions of these stars in the CMD are given in Fig. O 
and their properties are given in Table |2| We have also 
plotted the stars with the very best photometry, i.e., 3 581 
stars having the following quality-indices from ALLFRAME: 
\sharp\ < 0.25, |x - 1| < 0.25, ay < 0.04, and aj < 0.04. 
The sharp index measures how well the widths of the PSF 
and the observed star agree, x is essentially the standard 
goodness-of-fit index, and cry and cr/ are the estimated 
RMS errors on the PSF magnitudes. We refer to Sect. 4.1 
of Stetson et al. (2003) for a detailed discussion of the 
meanings of these indices. 

In Fig. 21 we present portions of the light curves for 
the ten stars with transit-like light-curve features in the 
NOTOl data set. The black points are V measurements 
while the grey points are / measurements. The times in- 
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Fig. 4. Light curves of transit candidates. Black and grey symbols are V and / measurements, respectively, and the 
noise level is given as cry and aj. The time to in each panel corresponds to the time HJD — 2 452 000 of the tick-mark 
at 0.0. Note that T6 has two transits-like events but of very different duration. 



dicated on the plots as to are HJD — 2 452 000 for the 
tick-mark corresponding to 0.0. 

We note that the telescope was offset slightly on the 
first night at the time to = 100.52 to the position which 
was used for the rest of the NOTOl campaign. This time 
coincides with the candidate transits seen in Tl, T2, 
and T3, and with the first of two candidate transits in 
TIO. Similarly, T6, T7, and T8 all have apparent transits 
around to = 106.60, but we note that the durations of 
these events are generally different. Finally, the observed 
transit-like event of T2, is very likely due to the star being 
close to a bad column on the CCD. 

We have only observed single transit-like events except 
for T6 and TIO. However, the durations of the two events 
seen for T6 are very different. T9 is found to be a long- 
period low-amplitude variable (V80, cf. Fig. [T)), so quite 
possibly the dip seen at to — 106.7 could be the beginning 
of an eclipse in a binary pair. This seems plausible from the 
position of T9 in the CMD in Fig.O The transit-like event 
seen for T4 is very odd, since the / measurements indicate 
a much longer transit duration than the V measurements. 

The most promising cases we have found are T5, T7, 
and T8, for which the transit events have a duration of 
1-3 hours. Furthermore, the depths are around lOmmag 
in both the V and / filters. These facts are not in contra- 
diction with the expected transit light curves of the short- 
period giant planets which are known from radial- velocity 
searches (Butler et al. I2002|l . 

As pointed out by Mallen-Ornelas et al. H2003|l . the 
shape of a transit caused by a planet will be characterized 
by steep ingress/egress slopes and a flat bottom during 
the transit. The precision of the present photometry is 



not good enough to sec the detailed shape of these shallow 
events. We would need further data of higher precision to 
confirm the transits and also to detect multiple transits. 

We finally note that the age of NGC 6791 is 8.0 ± 
0.5 Gyr according to Chaboyer et al. H1999|l (although we 
note that Stetson et al. H2003fl suggest the age could be 
significantly greater). Hence the expected rotational pe- 
riods of the stars at the turnoff would be of the order of 
weeks, and therefore star spots do not seem like a plausible 
explanation for the events we have found. 

To summarize, we cannot claim that what we have 
observed is indeed caused by a giant planet in any of the 
cases presented in Fig. 0| and Table O 

6. The Variable Stars in NGC 6791 

To search for variable stars we used the Stetson J-index 
(Stetson 1996J. To summarize, the J-index is a normal- 
ized sum of the deviation of each point (or pair of points 
with small separations in time) in the light curve com- 
pared to the expected noise level for that star. To deter- 
mine the expected noise level for a given star we used the 
mean repeatability of a large number of stars with similar 
magnitude. 

In the following we present the phased light curves of 
some of the variable stars in NGC 6791. In addition to 20 
known variables in our field we present 22 new variables. 
The locations in the CMD of the 22 new variables are 
shown in Fig.jS] The new variables all have low amplitudes 
and long periods. This explains why they were not found 
by Rucinski et al. ifTUMjl or Mochejska et al. iPTH^ due 
to the lower signal-to-noise ratios in their data sets. 
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Table 3. We give the name, ID number, a, 5, V, V — I, and B — I ioi the stars with transit-hke hght curves. All 
information except the adapted name in the first column is from Stetson et al. H2003|l . The positions of the stars in 
the CMD is shown in Fig. and the light curves can be seen in Fig. 0| As discussed in the text, the transits seen for 
Tl-4, T6, and T9 are not likely to be due to transiting giant planets (they are marked by a * symbol). 



Name 


ID 


Q2000.0 


feooo.o 


V 


V -I 


B-I 


Tl* 


4891 


19^20™ 43''.37 


+37°47'31"5 


17.675 


0.928 


1.818 


T2* 


6598 


19''20™-48':65 


+37°47'41'!l 


18.175 


0.979 


1.901 


T3* 


12616 


19''21™07':28 


+37°47'39'!9 


17.718 


0.915 


1.809 


T4* 


9049 


19''20'"55".47 


+37°44'06'!l 


17.493 


1.005 


1.839 


T5 


3671 


19''20™39':07 


+37°47'26''l 


17.509 


0.901 


1.784 


T6* 


3567 


19''20™38'!63 


+37°45'33''3 


17.635 


0.919 


1.811 


T7 


3723 


19*" 20'" 39". 25 


+37°45'39'!3 


18.002 


0.954 


1.827 


T8 


9020 


19''20'"55'!40 


+37°47'23''4 


17.029 


0.848 


1.670 


T9* 


12390 


19''21'"06':48 


+37°47'27'!8 


17.886 


1.015 


1.947 


TIO* 


8279 


19''20'"53':42 


+37°48'14"6 


18.907 


1.079 


2.092 




0.80 0,90 1,00 1,10 1.20 1,30 1.40 

V - Ic 



Fig. 5. Colour-magnitude diagram of NGC 6791 with 
standard V and V — Ic magnitudes from Stetson et al, 
H2()() 31, The 22 newly discovered variables have been 
marked. Only the stars with the best photometry are plot- 
ted to emphasize the location of the main sequence, the 
turnoff, and the red giant branch. 



the long period variables, our period estimates agree with 
previously published results. 

In Fig, IHI we show the phased light curves of some of 
the known variables and in Fig,|7|we present some the new 
variables. The ordinate gives the change in magnitude in 
the V filter and the phase is on the abscissa. 

The MO CHOI and SMR95 data are not plotted in 
Figs, El and Fig, [3 if the amplitude of the variable is com- 
parable to the noise level in these data sets. 

An interesting eclipsing binary that was previously 
known is V9, which was observed in all three data sets. 
It is classified as an RS Canum Venaticorum type star 
(RS CVn, see e.g.. Hall (UnEJ). The light curves in the 
V [left plot) and / filters are shown in Fig. |S1 where the 
"-}-" symbols are used for the NOTOl data, the filled grey 
symbols are the MO CHOI data, while the squares are the 
SMR95 data. The light curve of V9 is characterized by a 
deep eclipse which is modulated by a "distortion wave" 
probably due to dark spots on the surface of one of the 
stars in the binary pair. It is interesting to note that while 
the phase of the eclipse is unchanged the modulation wave 
has shifted phase from 1995 {SMR95) to 2002 {NOTOl and 
MOCHOl). The shape of the modulation wave seems to be 
unchanged, as was also noted by Mochejska et al. H2002fl . 

We finally note that the light curves presented here will 
be available from the Canadian Astronomy Data Centre 
(Stetson et al. lMT^ . 



In Table El we give a list of the variable stars we have 
observed. The columns in Table 01 contain the variable 
name, the ID number, right ascension, declination, period, 
the approximate amplitude in the V filter, and the variable 
type. The information in columns 2-4 is from Stetson et al. 

We have determined the period of each star by using 
the method of PERIOd98 fSDerl ll99^ . For some of the new 
long-period variables the time series is not long enough or 
the amplitude is too small to give a safe estimate of the 
period. In these cases we have only given an approximate 
period (indicated by the symbol). Except for some of 



7. Three Eclipsing Binary Stars in NGC 6791 

In Fig. El we show the light curves of three detached eclips- 
ing binaries. In each panel io is the time HJD — 2 452 000 
for the tick-mark corresponding to 0.0. The black and grey 
points are V and / measurements, while open circles are 
R measurements from MOCHOl. 

For each of the eclipsing binaries V20, V60, and V79 
we show the light curves from five nights, i.e., the nights on 
which the eclipses occur and one or more nights without 
eclipses for comparison. 

The eclipses of V20 and V60 were noted before by 
Rucinski l|1996(l and Mochejska et al. lj2002p . but now for 
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Table 4. The variables in the NOTOl field. The first column contains the variable name, where the new variables are 
printed with italics {V79-V100). Columns 2, 3, and 4 contain ID, right ascension, and declination from Stetson et al. 
H2003|l . Columns 5 and 6 are the periods and amplitudes (in the V filter) we have determined. The last column gives 
the variable star type. The positions of the variables in the colour-magnitude diagram are given in Figs. 01 and while 
the phased light curves of some of the stars are shown in Figs. El 13 andjHl 



Name 


ID 


0^2000.0 


^2000.0 


P [days] 


A [mag] 


Var. Type 


VI 


6271 


19''20'"47':61 


+37°44'32'!0 


0.2676766 


0.194 


EB, close 


V4 


8576 


19*" 20™ 541 22 


+37°48'23"8 


0.325667 


0.046 


EB, close 


V5 


5883 


19*" 20™ 46". 53 


+37°48'47"9 


0.312666 


0.017 


EB, close 


V6 


11376 


19''21'"02':72 


+37°48'49"1 


0.279062 


0.042 


EB, close 


V9 


6371 


19''20'"47':88 


+37°46'37'!4 


3.19174 


0.040" 


EB, close 


V16 


12695 


19''21'"07".59 


+37°48'09"6 


4.53 


0.026 


EB, close 


V31 


11307 


19''21™02".47 


+37°47'09"3 


3.34154 


0.015 


EB, close 


V33 


3886 


19''20™39'.81 


+37°43'54'^4 


2.35944 


0.040 


EB, close 


V38 


11652 


19''21™03".69 


+37°46'05"9 


3.845 


0.023 


EB, close 


V20 


8600 


19''20'"54':30 


-l-37°45'34"7 


14.4698 


0.3" 


EB, detached 


V60 


10746 


19''21'"00'!70 


-|-37°45'45"1 


7.510/7.1*' 


0.4/0.015'' 


EB, detached 


V79 


8943 


19''20™55".21 


-l-37°46'39"7 


~9.0038 


0.15" 


EB, detached 


V8(f 


12390 


19''21'"06".48 


-|-37°47'27"8 


4.447 


0.014 


EB? 


V14 


7637 


19''20'"51''.67 


-|-37°45'24"8 


5.4-5.8 


0.028 


Periodic Var. 


V17 


3626 


19''20'"38';88 


+37°49'04''6 


6.42566 


0.036 


Periodic Var. 


V41 


7397 


19''20™50'!97 


+37°48'24"8 


0.4807278 


0.039 


Periodic Var. 


V53 


10798 


19''21™00".84 


+37°44'35''4 


7.815 


0.039 


Periodic Var. 


V81 


3268 


19''20'" 49^.65 


+37°48'08"7 


7.58 


0.008 


Periodic Var. 


V82 


3857 


19*" 20™ 39". 71 


+37°47'36''l 


7.58886 


0.010 


Periodic Var. 


V83 


5848 


19''20™46".40 


+37°44'14"1 


7.014 


0.037 


Periodic Var. 


V84 


6305 


19*" 20™ 47". 71 


+37°44'58"2 


1.6333 


0.03 


Periodic Var. 


V85 


7011 


19^20™ 49".86 


+37°45'50"9 


4.100 


0.007 


Periodic Var., EB? 


V86 


7099 


19*" 20™ 50". 13 


+37°48'31"7 


7.39 


0.011 


Periodic Var. 


V87 


8033 


19*" 20™ 52". 78 


+37°44'58"8 


6.971 


0.006 


Periodic Var. 


V88 


8088 


19''20™52".91 


+37°46'36"9 


6.91 


0.008 


Periodic Var. 


V89 


9432 


19*^20™ 56".64 


+37°46'36"2 


5.7617 


0.028 


Periodic Var. 


V90 


10169 


19''20™58".86 


-|-37°44'47"2 


5.71 


0.008 


Periodic Var. 


V91 


10699 


19''21™00"54 


+37°48'40"7 


5.13 


0.011 


Periodic Var. 


V92 


11380 


19''21™02".73 


+37°46'00"8 


7.24 


0.008 


Periodic Var. 


V93 


12049 


19''21™05".23 


+37°47'08"5 


0.94928 


0.007 


Periodic Var. 


B7 


12652 


19''21™07".40 


+37°47'56''5 


~13.3 


0.15 


Long period 


V62 


11475 


19''21™03".06 


+37°43'51"8 


~17 


0.027 


Long period 


V65 


7920 


19''20™52".47 


+37°47'30"5 


8.613 


0.008 


Long period 


V67 


11645 


19''21™03".67 


+37°48'03"7 


~10 


0.05 


Long period 


V73 


8108 


19''20™52".97 


+37°46'52"5 


~20.9 


0.02 


Long Period 


V94 


4640 


^gh 20™ 42". 50 


+37°44'36"9 


~15.57 


0.01 


Long period 


V95 


4805 


19*^20™ 43".05 


+37°47'32"5 


9.477 


0.012 


Long period 


V96 


5485 


19*" 20™ 45". 26 


+37°45'48''8 


9.307 


0.012 


Long Period 


V97 


6777 


19^20™ 49".17 


+37°49'14"8 


9.563 


0.013 


Long period 


V98 


9353 


19''20™56".42 


+37°45'38"4 


~8.99 


0.01 


Long period 


V99 


9568 


19''20™57".06 


+37°48'12"2 


~10 


0.003 


Long period 


VI 00 


11111 


19''21™01".81 


+37°45'42"0 


12.5215 


0.01 


Long period 



(a) For the eclipsing binaries V9, V20, and V79 the depth of the primary eclipse in the V filter is given. 

(b) V60 is variable with a low amplitude and a period of P ~ 7.1 days while the period estimated from the 
secondary and primary eclipses is 7.510 ±0.005 days. The primary eclipse depth is 0.4 mag in V. 

(c) V80 was detected in our search for transit-like events and is also given the name T9 (cf. Fig. 0J. 
Furthermore, V80 is variable with a low amplitude. 



the first time we have enough time coverage to be able to 
constrain the periods. From the difference in time between 
the eclipses we find periods of Pbin — 14.470 ± 0.001 and 
7.510 ±0.005 days for V20 and V60, respectively We note 



that V60 is variable with an amplitude of 0.015 mag in V 
at a period of Pvar — 7.1 ± 0.2 days. 
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Fig. 6. Phased light curves in the V filter for some of the known variables in NGC 6791. In most panels only the 
NOTOl data (+ symbols) are plotted due to the much lower noise level compared to MOCHOl (triangle symbols) and 
SMR95 (box symbols). Notice the change in the light curve of V4. 



7.1. The Detached Binary V20 

We have been granted observing time to obtain more de- 
tailed photometric coverage of the primary and secondary 
eclipses of V20, and we plan to carry out spectroscopic 
measurements to obtain the radial velocity curve. With 
these observations it will be possible to constrain the 
masses, radii, and luminosities of the binary components 
of this system. 

As will be demonstrated in Sect. 17.21 the primary com- 
ponent of V20 is very close to the turnoff. From the know- 
ledge of the mass of the primary component, we expect 
to be able to determine the age of NGC 6791 with good 
accuracy. This can be done from the relation between age 
and mass of the turnoff stars which will rely on theoretical 
isochrones. 

This has recently been done for the globular cluster 
w Gen by Thompson et al. (20011 and Kaluzny et al. 
for eclipsing stars around the turnoff. Their masses 



could be determined with an accuracy of just 1 %, includ- 
ing systematic errors. In this case one of the systems is 
right after the turnoff and Kaluzny et al. (|2002|l were able 
to constrain the position of the theoretical isochrones and 
in turn estimate the age of uj Cen to within 0.6 Gyr — at 
least a factor of two better than any previous age deter- 
mination for a globular cluster. 

7.2. Estimating the Parameters of the Detached 
Binary V20 

In Fig. El the colour- magnitude diagram of NGG 6791 is 
shown. The location of the eclipsing binary V20 is shown 
with a star symbol above the turnoff region. In the fol- 
lowing we will describe how we have determined the pa- 
rameters of the two components that make up the binary 
V20. 

We have fitted a fiducial locus to the turnoff region 
shown as the solid "primary" curve in Fig.^l and in the 
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Fig. 7. Phased light curves of some of the new low-amplitude variables. Only the NOTOl data are plotted. The ID 
number is given in the top left corner of each panel and correspond to Stetson et al. H2003|l . The period (in days) is 
given in the lower left corner. More details are given in Table 0] 



V09: 






V09: 






0.1 










"^""^ 




% 00*° 


0.0 


2 ilPrft 




□ - 


























0.1 




K 






3.192d 


tj - 

4° = 







Fig. 8. 

black 4 
(filled g 
in / 



,0 0.4 0.8 1.2 0.0 0.4 0.8 1.2 

Phased light of the RS CV-type variable V9 in NGC 6791 in the V {left plot) and / filter {right plot). The 
symbols are data from NOTOl while the grey symbols are from either SMR95 (box symbol) or MOCHOl 
xey symbols) The MOCHOl data are not plotted in the right plot since Mochejska et al. H2002|l did not observe 



12 



H. Bruntt et al.: A Search for Planets in the Old Open Cluster NGC 6791 



0,0 
-0.1 

ai 
a 

E -0.2 
<: 

-0.3 

-O.A 














< 

t.=-aoa9.60 


t„=-yOBy.45 


to=10y.55 


* • o 

w 

t.= 109.80 


V20 

t.= 106.60 



-0.1 0.1 
At 



G G 
-0.2 

a 
E 
< 


> 


t 






V60 


-0.4 












-0.6 


t„=-aoa8.40 


t„=10y.80 


to=104.50 


to=105.50 


t„=106.50 




-0.1 0.1 
At 


-0.1 0.1 


-0-1 0-1 


-0-1 0-1 


-0.1 0.1 



0-05 

00 .v*>iir>' 



CP -0-05 

I -0-10 

-0-15 
-0-20 
-0-25 . 



t„=-20t3a-40 



t„=-2083.40 : t„=104.40 to=105.50 to=106.50 



V79 



Fig. 9. Light curves of three edipsing binary stars. The 
time corresponding to 0.0 on the abscissa is HJD — 
2 452 000 and is given at the bottom of each panel as ^o- 
One or more nights without eclipses are shown for compar- 
ison. Black and grey symbols are V and / measurements, 
while the open circles are R measurements from MOCHOl. 



" 1 f 



19 



20 




1 .0 



1,2 1.4 
V - I 



1 .5 



Fig. 10. The I vs. V ~ I colour-magnitude diagram of 
NGC 6791 on the standard system. The location of the 
eclipsing binary V20 above the turnoff is shown with a star 
symbol. To find the possible positions of the two compo- 
nents of the binary V20, we assume that the primary star 
is on the solid black curve and calculate the corresponding 
colour and magnitude for the secondary component which 
is given as the dashed line. For three points we have shown 
the location of the primary and secondary star, i.e., for 
the primary star being at the turnoff (TO), further down 
the main sequence (MS), and the point where the dashed 
curve intersects the main sequence (Bin). 



following we assume the primary component will be on 
this curve. For each point on the "primary" curve we cal- 
culate what the magnitude and colour of the secondary 
component should be, while requiring that the combined 
colour and luminosity of the primary and secondary com- 
ponent must yield the observed position of V20. The in- 
ferred position of the secondary component is shown as 
the dashed "secondary" curve. 

For three positions on the "primary" curve we have 
marked the corresponding position on the "secondary" 
curve: at the turnoff (TO), on the main sequence (MS), 
and the point where the main sequence and the "sec- 
ondary" curve intersect (Bin). The latter is the best fit of 
the location of the two components making up the binary 
V20. Assuming that this inference is correct, their colours 
are determined to within less than ±0.01 mag; these are 
given in Table We will now use these colour estimates 
to infer the properties of the components of V20. 

We have used the calibration of Tog vs. {V — I)j 
for main sequence stars with solar metallicity by Alonso 
et al. I|1996|) . The calibration of Alonso et al. (1996) re- 
fer to Johnson Ij magnitudes, and so we have used the 
calibration of Fernie H1983fl to convert from the Kron- 
Cousins system we have used, i.e., {V — T)j ~ —0.005 + 
1.273(1^ — I)c. For the interstellar reddening we use 
E{B -V) = 0.105 ± 0.025 (Chaboyer et al. ITTMIl and 



using Taylor l|TM)|l we find E{V - I) = 0.13 ± 0.03. We 
note that the error given on E{B — V) does not include 
systematic errors which could very likely be of the order 
a{E{B — V)) ~ 0.05. However, the analysis of the colour- 
magnitude diagram of the cluster by Stetson et al. H2003|l 
suggests a reddening E{B — V)'^ 0.09, entirely consistent 
with the Chaboyer value that we employ. 

From these numbers we can determine the de-reddened 
magnitudes and colours and thus estimate the Tgff's of 
the two components of V20 to within 200 K, includ- 
ing the systematic uncertainty in E{B — V). The re- 
sults are given in Table El We obtain the bolometric cor- 
rection from Bessell et al. p998|l for atmospheric mod- 
els with solar metallicity and no overshooting. We can 
then determine the relative luminosities of the stars using 
Mboi;p-Mboi;s = {Vp + BCp)-{V, + BC,) = -2.84 ±0.14 
which is equal to —2.5 log(Lp/Ls), and hence the luminos- 
ity ratio is ip/Ls = 13.7 ± 1.8. When using L oc R^T^s^ 
we find a radius ratio of R^/Rp = 0.41 ± 0.08. 

We have also estimated the relative luminosities of the 
binary components by modelling the observed light curve 
of V20 following the method outlined in Clausen et al. 
(2003). We assumed Tcs = 5 200 K for the primary com- 
ponent and from the modelling of the light curves using 
only the light curve in V we find Tcff = 4 400 ± 250 K for 
the secondary component. The large error in T^g is mainly 
due to the depth of the secondary eclipse being different 
by ~ 0.02 mag in the SMR95 and NOTOl data set. From 
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the light curve modelling we also derive the radius ratio 
to be Rs/Rp = 0.55 ± 0.1. Both T^g and the radius ratio 
are in good agreement with what we have found by the 
two methods described here. 

Obviously, much more detailed photometric coverage 
of the primary and secondary eclipses of V20 are needed 
to further constrain the properties of the two components. 
As already mentioned, we need spectroscopic observations 
of the complete orbit to be able to constrain the masses 
of the components. 

8. Conclusion and Future Prospects 

We have described the search for planet transits in the 
open cluster NGC 6791 from seven nights of observations 
with the Nordic Optical Telescope. The photometric preci- 
sion is adequate for detecting transits by close-in Jupiter- 
sized planets. We have found ten stars with transit-like 
events in the light curves, but we argue that most of these 
are probably not caused by a transiting giant planet. The 
observed transit events may be due to instrumental effects 
(bad columns, offset of the telescope) and in one case there 
is evidence that the star (T9) belongs to a binary system. 
For the remaining three transit-like events the precision 
(3 mmag) is not good enough compared to the very shal- 
low transit depth (~ 10 mmag) to tell whether the transits 
are really caused by giant planets, e.g., to sec whether the 
light curves are flat during the transit. To confirm the 
validity of these three planetary transit candidates the 
observation of multiple transits with higher precision are 
required. 

In July 2002 Piotto, Stetson et al. have carried out 
a more ambitious multi-site campaign on NGC 6791 to 
search for planet transits. We have obtained observations 
on eight at least partially clear nights of a ten-night run 
on the CFHT 3.6 m telescope in Hawaii during 2002 July 2- 
12. In the same period of time observations from Mexico 
and Italy were also carried out. The weather was not opti- 
mal during the multi-site campaign on NGC 6791 in 2002, 
but we have a higher signal-to-noise ratio compared to the 
NOTOl data and also a higher duty cycle over 10 nights. 
This will improve the statistics for detecting multiple tran- 
sits. The reduction of data from the 2002 campaign data 
is now under way. 

In addition to our own data from NOTOl we have used 
data from two previous photometric campaigns. The pho- 
tometric precision of these other data sets is not high 
enough to confirm the shallow transits we have found. 
We have used the complete data set to redetermine the 
periods of the 20 known variables in the NOTOl field. 
Furthermore, we have discovered 22 new long-period stars 
with low amplitudes. 

We have determined the periods of two detached 
eclipsing binaries, V20 and V60. From the position of V20 
in the colour-magnitude diagram we are able to constrain 
the temperatures of the binary components. This is in 
good agreement with results from the modelling of the 
primary and secondary eclipses. We find that the primary 



component is located just after the turnoff, while the sec- 
ondary component is a much fainter main sequence star. 
We plan to make thorough photometric and spectroscopic 
observations of V20 to constrain the masses of the binary 
components. From the mass of the primary star we can 
determine the age of NGC 6791 from a comparison with 
the calculated turnoff mass for isochrones of different ages. 
This will give an age estimate independent of interstellar 
reddening. 
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Table 5. The derived parameters for the two components of the binary star V20. We found the colours from Fig. IIUI 
and in the text we explain how Toff was determined. The estimated errors are given in parenthesis. 
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